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We study a simple transport model driven out of equilibrium by reservoirs at the boundaries, 
corresponding to the hydrodynamic limit of the Symmetric Simple Exclusion Process (SSEP). We 
show that a non-local transformation of densities and currents maps the large deviations of the 
model into those of an open, isolated chain satisfying detailed balance, where rare fluctuations are 
the time-reversals of relaxations. We argue that the existence of such a mapping is the immediate 
reason why it is possible for this model to obtain an explicit solution for the large-deviation function 
of densities through elementary changes of variables. This approach can be generalized to the other 
models previously treated with the Macroscopic Fluctuation Theory 



Systems driven by reservoirs at the borders are possi- 
bly the first step in complication as one enters the non- 
equilibrium realm. They are already not solvable in gen- 
eral, but many elegant and striking results have been 
obtained in recent years (for a review, see Q and refer- 
ences therein). Recently, Derrida, Lebowitz and Speer Q 
(DLS) obtained an exact expression for the large devia- 
tion function of the density profile in the Symmetric Sim- 
ple Exclusion Model (SSEP) - a one-dimensional trans- 
port of particles - using a matrix method that had been 
developed previously [4j|. This is a major achievement, as 
large deviation functions are a natural extension of free 
energy to out of equilibrium situations. 

For the sake of generality, it is desirable to have a 
purely macroscopic approach, that does not rely on the 
strong symmetries of the microscopic model. This was 
done within the Macroscopic Fluctuation Theory 0, [|| . 
This is a version of WKB theory, valid in the limit of large 
coarse-graining scale N: in the usual manner of Semi- 
classical Theory and Geometric Optics, the logarithm of 
wavefunctions evolving with the Fokker Planck operator 
obey a Hamilton-Jacobi equation whose characteristics 
are trajectories satisfying Hamilton's equations. 

For the hydrodynamic limit of the SSEP, Bertini et 
al.Q were able to integrate explicitly the correspond- 
ing Hamilton-Jacobi equations and recovered the large- 
deviation function. They thus showed that such a so- 
lution is in principle logically independent of the exact 
microscopic solution. Their derivation amounts to a se- 
ries of carefully chosen changes of variables, each of one 
yielding a miracle, whose cause and degree of generality 
one may still wish to ascertain. 

In this paper we uncover a transformation of the large 
deviations of densities and currents of the original driven 
model into those of a 'dual' isolated, equilibrium chain. 
Large deviations and optimal trajectories are easily ob- 
tained in this representation using the detailed balance 
property, and can then be mapped back to the original 
setting. Moreover, emergences of rare fluctuations in the 



dual model are the time reverse of relaxations to the av- 
erage profile, but this symmetry is lost in the mapping 
back to the original model. This accounts in this case for 
the lack of Onsager-Machlup symmetry [6| between birth 
and death of a fluctuation, which has received consider- 
able interest [H, 0, 0] over the past few years. 

In what follows, we shall first rephrase the Macroscopic 
Fluctuation Theory in a path integral language, then sin- 
gle out the specific properties which enable one to com- 
pute the large deviation function for the system ([I]) and 
finally encapsulate them in the existence of a dual model 
with detailed balance. 

Hydrodynamic Model 

We consider the fluctuating hydrodynamic limit of the 
SSEP in contact with two reservoirs [9|] 

P = -VJ; J = ~Vp- y/a^T) ^ 

p(0,t)=p o ; p(L,t) = p L 

Here 77 is a white noise of variance 1 /N and a = p(l—p). 
The probability of a trajectory is given by [lfj 

P ~ Jv[p lV }5(p + VJ)e- N f dtdx £ 

~ J T>[p, p,r)]e~ N 1 ' dtdx {p{f>-T! A P-v(V^v)} + 4} 

After an integration by parts (which shall be frequent 
and unannounced in the following), the integration over 
X] gives 

Jv\p,p]e- NS ^ = J V[p,p]e- N i dtAx ^P- n ^ (2) 

where we have introduced the Hamiltonian density de- 
fined by H = j'Hdx: 

U=\ [cT p Vp 2 + pAp] (3) 
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Large Deviations: General strategy 



The first shift maps 



In the large- N limit, the probability of observing a profile 
p*(x) in the steady state scales as P(p*) ~ e^ NF ^ p *^. 
We wish to calculate the large-deviation function T ', 
which is given by the action of the 'instanton'. This is 
the trajectory starting in a neighbourhood of the station- 
ary profile p that extremises the action ([2]). converges to 
p* (x) at a large time T (and thus has at all times H = 0) 
and satisfies the spatial boundary conditions. The prob- 
lem then reduces to solving the equations of motion: 



P = 



SH 
5p(x) 
SH 
Sp(x) 



1 



Ap-V[a p Vp] 
-^+(2,-1)^ 



(4) 



with the space and time constraints 
p(x,0) = p(x) = po + x- 



p(x,T) = p*(x) 



L ' r ^' (5) 

p(0,t)= Po ; p(L,t)=p L ; p(0, t) = p(L, t) = 

The last equality simply says that no fluctuations are 
allowed at the contact with the reservoir [l], 0]. Alter- 
natively, one can solve the classical equations ((4]) via the 
Hamilton- Jacobi formalism [13j | and that amounts to the 
strategy followed by Bertini et al. 

Mapping the problem into a downhill one. 
Detailed balance 

The fact that equations ((4]) derive from a stochastic prob- 
lem implies that there is a family of explicit 'downhill' 
(zero noise) solutions: p(t) = C st and p = |Ap. For 
C st = 0, the corresponding action S[p, p] is zero. Such a 
solution is not what we are looking for, as it relaxes into 
and not out of the stationary state and hence does not 
satisfy the boundary conditions in time ([5]). 

A strategy to find solutions of equations like (|4]) is to 
make a change of variables that maps the original prob- 
lem into another one of the same form - i.e. that for- 
mally derives from some other stochastic problem, but 
such that the downhill solutions of the new problem obey 
the correct boundary conditions in time and in space [1 51 ] . 

For a chain at equilibrium, a simple procedure can be 
followed to do so, taking advantage of the detail balance 
relation of Hamiltonian ([3|) , which can be made explicit 
by writing 



H = \vp a p 



log- 



P 



\^P ° P V 



SVp 
Sp 



with V p = J ddplogp + (1 - p) log(l - p)\. At the level 
of the action [l6| , detailed balance means that this form 
is left invariant by a succession of two transformations 



p -> p + log 



1-P 



SVp 
6p 



(p,t)^(-p,T-t) (6) 



n^n = -vp u p v 



sv 



S[p, p] 



w, 



/'JO 



Sp 

dtdx | p'p — H | 



(7) 



The new H has the form of a stochastic problem, and the 
equations of motion associated to ([7]) admit a 'downhill' 
solution p — C st which in the old variables reads 



P = - 2 A P 



P = log 



1-P 



(8) 



and corresponds to the optimal uphill trajectory. The 
second shift in ^ shows that this trajectory is the time 
reverse of a diffusive trajectory. For a chain driven out 
of equilibrium by the boundaries, this simple strategy 
fails: (jSJ) is not compatible with the spatial boundary 
conditions ([5]) for p. 

Before going on, let us note a striking property, specific 
to the problem |T]). Rearranging the Hamiltonian density 



n=-p 



(Vp)' 



P-i- 



(Vp) 2 



= ~P- 



(Vpf 



p-l-f 



(9) 



where Vp — J dx[log Vp], Ti. can formally now be seen as 
deriving from another stochastic dynamics: p = — ^Ap — 
^(Vp) 2 +rjVp with a further detailed balance symmetry, 
induced by the two transformations 

P^P+^; (p,t)^(i-p,T-t) (10) 
Sp 

As in the previous paragraph, two further classes of so- 
lutions can be directly read in |9]): p = 0, p = — ^Ap — 
\Vp 2 and p — 1 = p = iAp + |Vp 2 , respectively. 
However, none of these solutions satisfies the boundary 
conditions. This additional symmetry is a signature of 
the existence of the dual model, as we shall see below. 
This is the first time the specific form of the model (in 
particular its one dimensional nature) plays a role. 

A Solution for the classical problem 

Let us now paraphrase Bertini et al., following a con- 
cise but unrigorous Hamiltonian - rather than Hamilton- 
Jacobi - approach. To make contact with the exact so- 
lution, we rewrite the intermediate variable p appearing 
in equation ([7]) in terms of the DLS variable F defined 
by F = (1 + eP)^ 1 . To keep the action in a Hamilto- 
nian form, we also introduce the canonically conjugate 
variable F = ^ . This maps the action into 



F(l-F) 



S = dx 



p log p + ( 1 - p) log(l - p) + p log 



1 



1 - F 



F 



-!- / d.rdl { FF + ^FVF 2 







\*- 2 " 




1 -F_ 
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This expression can be further simplified by making the 
shift F -> F + to get 



S 



plog£ + (l_p) log^— J; 



Ax 

dxdt { FF 



1 F 2 (VF) 2 + ^VFVF 



(11) 



2 v ' 2 
The overall mapping from the initial (p, p) to (F, F) reads 



F : 
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P +(l-p)eP- 



F=(l-p)(e/-l)-p(e-P-l) (12) 



The boundary conditions J5]) are now given by F(0, t) = 
po, F(L,t) = p L and F(0,t) = F(L,t) = and the 
equation of motions are 



F = ^AF-F(WF) 2 ; 



F = — -AF - V 
2 



F 2 VF 



(13) 



Modulo an integration by part, the last integral in fill) 
is of the form J dxdt |fF — Hp| where 



H 



-F 



(VF) ; 



F 



AF 

WW 



-F 



(VF)' 



F 



SV F 
5F(x) 



This is very similar to ([9]), but with Vjr = J dx ln(VF). 
Quite surprisingly, we have once again obtained an action 
formally deriving from a stochastic dynamics, satisfying 
the detailed balance symmetry induced by 



F -> F 



AF 
WW 



(F,t) -> (-F,T-t) 



(14) 



As in all the previous examples, two classes of solutions 
are immediately available. First, F = 0, F = ^AF cor- 
responds to a downhill diffusive solution, equivalent to 
p = ^Ap and which is not what we are looking for. 
Instead, an uphill trajectory is provided by requiring 
F = ■ Together with the equation of motion (fT3")) . 
it implies 



F 



:AF: 



F + F(l -F) 



AF 
VF 1 



(15) 



Amazingly, this time, its solution satisfies both spatial 
and temporal boundary conditions, as is easy to check. 
The corresponding action is the large deviation function 



T = dx 



plogj; + (l-p)10g^— ■ 



Dual Model 



1 T 



log VF 



(16) 



The above derivation consists of changes of variables that 
read like a sequence of miracles, at the end of which 
we are able to find an 'uphill' solution that satisfies the 
boundary conditions. In the following we show that all 
the surprises can be though of as deriving from only one: 



the action for a chain in contact with two reservoirs can 
be mapped, at the level of large deviations, to that of an 
open, isolated chain. Starting from the action (fTT|) . we 
introduce the non-local variables 



F' = VF; F = V 
which takes the action into 
S= [dx 



F>-1 
F' 



VF' 



VF' 

- 2 

F' 



(17) 



1 T 



plog | + (1 - p) log |-J; + logF' - F'F' 



dxdt \ F'F' + iF' 2 (VF') 2 + ^VF'VF' 



Remarkably, the form of this action is, up to boundary 
terms, the same as (fTTj) . This suggests that we complete 
the mapping 



(P>P) 



fp,F) no ^i°, cal (F',F' 



(p',p' 



(18) 



where the relation between (p',p') and (F',F') is of the 
same form as (IT2l. One then obtains 



S= dx 



p log | + (1 - p) log — J; + log VF 



^-p'log^-(l-pOlogl^ 
I — F' F' 1 — F' 



(19) 



S' 



where S' = j dtdx {p'p' — H'} and H' is formally equiv- 
alent to @ but for the primed variables: 



1 



1, 



5"= / dtdx jp'p' - -oy Vy + -Vp'Vp' 



2 " - r '2 
The overall change of variable (fT5|) , which reads 
1 



(20) 



V 



1-eP 
P 



p+(l-p)eP 



1 - p{e p + e~ p - 2) 



l-p'(e p ' + e ~P' -2) 



(21) 



thus maps the action of the hydrodynamics limit of the 
SSEP into another SSEP. We shall next show that the 
boundary conditions transform in such a way that the 
dual chain is isolated. Before going on, it is instructive 
to introduce the classical spins S' 

S' z = 2p' - 1 S' + = 2(1 - p')eP SL = 2p'e-P (22) 

This is the usual connection between spin chains and par- 
ticle models [HI, [HI in the hydrodynamic limit. The 
Hamiltonian density reads H' = -|VS' • VS' . It is 
invariant by simultaneous rotation of all the spins, which 
means that to the three 'charges': 

2p'-l = S^ ; F' = S' z +iS' ; F'(1-2F') = S' x -1 
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correspond three currents 
J P > =~Vp' + <r P >Vp' 5 Jp, = \^F' + F 2 VF' 
Jp { i-2F>) =(! - 2F')^f + [F'+F'(l 2F')}VF> 



(23) 



that are conserved in the bulk. The boundary conditions 
read for the primed variables 



F>-1 
F' 



= 0; 



x=0,L 



F' = pi — po (24) 



Note that all trajectories satisfying both the equations 
of motion (|13p and the boundary conditions satisfy also 
AF = V-F' = on the boundaries. Together with the 
l.h.s. of (|24| and the definitions (fT2| . this implies that 
VF'jVF'jVp' and V/5' vanish separately at the bound- 
aries, as do, consequently, all the currents (f2"3"| : the dual 
model in the primed variables is an isolated chain. This 
condition alone, supplemented with the r.h.s. of ([24]) . 
encompasses all the original boundary conditions. 

Let us now turn to the trajectories. The stationary 
profile p maps to a flat profile p' , but whereas F' is 
constrained by the r.h.s of (f2"4"|) , the precise value of p' 
is arbitrary. Relaxations: the diffusive trajectories of 
the initial model correspond to p = \Ap, p = 0, and 
thus F — 0, as can be read in (fT2"|) . Together with 
([2"Tj) . this implies V/5' = 0. As the primed variables 
also satisfy the equations of motion (H|), the resulting 
trajectories evolve with p! = \Ap' . Relaxations thus 
map into relaxations. Excursions: The instanton equa- 
tions (HSl) imply VF' = (cf. ([TT|)). Using the relation 

(1 1 21) f^iT" fVlo T-»yimo/-l ri a VJ o fViic Triolrlc V7 rj — ^ P — 



for the primed variable, this yields V p' — = 0. 
Together with the equations of motion ((4]), this shows 
that the densities evolve with p' = — ^Ap': excursions 
are the time reverse of relaxations in the primed vari- 
ables. The action S'[p' , p'\ of such an uphill trajectory is 
J &x[p' logp' + (f - p') log(l - p')]q . As F' is constant 

along the instanton and p' is a constant of motion, the 
overall action (fH?)) reduces to the large deviation function 
(HHJ) as it should. 

Let us now show how the existence of the primed model 
accounts for the many non-local symmetries of the ini- 
tial one. The action of the dual, isolated, chain is invari- 
ant under a composition of any simultaneous rotation 
and reflection of all the spins, followed by time-reversal. 
Written in the original variables, this corresponds to a 
group of spatially non-local symmetries of the out-of- 
equilibrium chain. Transformations of the primed model 
which do not conserve F' = PL ~ Po correspond to 
symmetries of the initial model broken by the bound- 
ary conditions. This was for instance the case of (fl~0|) . 
If we further impose the conservation of J Q F, we are 
left with a transformation that in spin representation 
amounts to (S' x , S' y , S' z , t) -> (S' x , S' y , S' z , T - t). This 



is in the primed variables the non-local mapping (fT4| be- 
tween 'downhill' diffusive solutions and the instantons of 
the initial problem. 

Conclusion 

In this letter we have shown that the remarkable proper- 
ties of the hydrodynamic model allowing for its direct, ex- 
plicit solution are attributable to the existence of a dual, 
equilibrium model. Whereas the derivation we have pre- 
sented is very specific of the model we studied, it can eas- 
ily be extended to the other cases recently solved within 
the Macroscopic fluctuation theory, including the Kip- 
nis Marchioro Presutti model A further extension 
to a larger class of non-equilibrium systems or to micro- 
scopic models is still an open and challenging question. 
The symmetry between excursions and relaxations in the 
dual model is broken by the mixing of variables due to the 
non-local mapping (| 1 T|) . ft would thus be very interesting 
to see if one can construct non-local quantities (like the 
currents of the dual model) which would be symmetric, 
and measure them experimentally. 
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Conversely, if one knows the large deviation function 
J-(p) , such a transformation is given by p —* p + ^£ [l|] . 
In this paper we shall only invoke detailed balance at the 
large iV 'classical' level 

The equations of motion read S = ^AS A S. 



